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Abstract In this work we quantify the degree to which star-forming events are 
clumped. We apply a precise and accurate technique to calculate the correlation 
dimension of the distribution of HII regions in a sample of disk galaxies. Our 
rehable results are distributed in the range 1 .5 <; Dc <; 2.0. We get significant varia- 
tions in the fractal dimension among galaxies, contrary to a universal picture some- 
times claimed in literature. The faintest galaxies tend to distribute their HII regions 
in more clustered (less uniform) patterns. Moreover, the fractal dimension for the 
brightest HII regions within the same galaxy seems to be smaller than for the faintest 
ones suggesting some kind of evolutionary effect. 



1 Introduction 

There is clear evidence that gas observed in external galaxies (both irregulars and 
spirals) follows hierarchical and self-similar patterns 1201 [121 1211 [3l |5]| . This frac- 
tal picture is consistent with the distribution of star fields and star-forming sites on 
galaxy-wide scales Q [15] |6l IH [H 121 . However, it is not clear whether this kind 
of fractal distributions are connected/related or not to some properties of the host 
galaxies, such as radius, rotation, brightness, morphology, etc. In spite of the great 
variety of fractal dimension values reported in the literature for different galaxies 
(for both the gas and the distribution of star forming sites) most of the authors argue 
in favor of a more or less universal picture. In this universal description, the con- 
stancy of the fractal dimension is a natural consequence of the fact that the same 
physical processes are structuring these systems. However, there are some indica- 
tions that the situation could be more complicated. The fractal dimension of the 
distribution of HII regions could be different in grand design and flocculent galax- 
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ies IfTTI . and the brightest galaxies could have fractal dimensions higher than faintest 
ones ifTSl [131 . On the contrary, Q do not find any correlation between the fractal 
dimension and the galactic properties, but their uncertainties are so large that the 
robustness of this conclusion is questionable. 

Part of the problem that prevents achieving unequivocal conclusions lies in the 
great diversity of analysis techniques used in the literature and/or the application 
to not large enough samples of galaxies. Here we use a carefully designed method 
that has been tested previously on simulated data and that clearly establishs its ac- 
curacy and applicability depending on the sample itself. We apply this method to 
the most complete sample of galaxies that we have found in literature expecting to 
draw significant conclusions regarding this matter. 



2 Fractal dimension calculation 

We use the correlation dimension fTOl which gives robust results when dealing with 
distributions of points in space. First we have to calculate correlation integral C(r) 
which represents the probability of finding a point within a sphere of radius r cen- 
tered on another point. For a fractal set C{r) scales at small r as 

C(r)^A , (1) 

being Dc the correlation dimension. When evaluating C{r) for real data this power- 
law behavior is valid only within a hmited range of r values. For relatively small or 
high r values C{r) deviates from the expected behavior and Dc tends to be underes- 
timated. We have developed an algorithm (see details in ifTSl and |[T9l ) to calculate 
Dc which implements objective and suitable criteria to avoid both boundary effects 
and finite-data problems at small scales. Moreover, the algorithm estimates an un- 
certainty associated to Dc by using bootstrap techniques {Oboot)- 

An important point is the dependence of the measured dimension on the sample 
size because many times the number of available data points is rather small. Figure[T] 
shows what happens when the calculation is done on the same set of simulated 
fractals but randomly removing points to reach smaller sample sizes A^. Since we are 
dealing exactly with the same original data, the observed decreasing in the average 
Dc value with decreasing have to be attributed exclusively to sample size itself. 
Thus, our algorithm is able to estimate the fractal dimension in a reliable way when 
applied to random subsamples of fractal distributions of points, but if the sample size 
is too small (A^ 100) a bias tending to underestimate the dimension is produced. 
The existence of this bias is related to the random variations in the simulated (or 
observed) data set and not related to the uncertainties in the determination of Dc 
iPboot)- The latter ones are always smaller than the error bars in Figure[T] 
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Fig. 2 Left panel: positions (units are in kpc) of the HII regions in NGC 6946. The dashed line 
circle indicates the radius R25. Right panel: the calculated correlation integral C(f ) for NGC 6946. 
The lower limit of the best linear fit (solid line) is indicated by the vertical aiTow. 



3 Results 



We have used VizieR and ADS databases, in conjunction with the papers IS) and 
H, to search for catalogs of external galaxies containing positions of HII regions 
available in machine readable format. We found a total of 93 spiral galaxies with po- 
sitions for at least 50 HII regions. We have also included data for 8 irregular galaxies 
with HII positions provided by D. Hunter ifTTll . The properties of the selected galax- 
ies are listed in Table 2 of |fT9l available through VizieR at CDS. This table contains 
the galaxy name, Hubble type, position ^ and inclination / angles, morphological de 
Vaucouleurs type T , arm class A, distance D, B-band absolute magnitude Mb, radius 
of the isophotal 25 mag/arc sec^ R25, and the maximum rotation velocity Vrot- The 
deprojected positions of the HII regions for the sample of 101 disk galaxies, which 
can be seen in Figure 3 of |fT9l , are also available through VizieR. 

Figure|2]shows the result for one example galaxy. The left side panel in this figure 
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Fig. 3 Calculated dimension 
D^. as a function of the number 
of available data A'(HII) 
for spiral galaxies (solid 
circles) and irregular galaxies 
(open circles) in the sample. 
The error bars indicate the 
uncertainties obtained from 
bootstrapping. 
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shows the distribution of HII regions in NGC 6946. The right side panel shows the 
corresponding correlation integral C{r). The vertical arrow indicates the value of 
r below which C{r) is poorly estimated ifTsl [T9l . The rest of the data exhibits a 
characteristic fractal behavior The slope of the best linear fit gives Dc, in this case 
L», = 1.75. 

The result of calculating the correlation dimension of the distribution of HII re- 
gions in the full sample of galaxies is shown in Table 3 of lfT9l (also available via 
VizieR). Figure [3] shows the calculated correlation dimension for all the galaxies 
in the sample as a function of the number of available data points. We see that as 
A^(HII) decreases, the uncertainties increase and the obtained dimensions are sub- 
stantially more spread out toward lower values. This trend is very similar to the 
behavior shown in Figure [T] It is clear that at least part of this trend is due to a bias 
in the estimated value of Dc for galaxies with small number of HII regions. To over- 
come this bias, we focus our analysis on galaxies having more than 200 HII regions 
(46 spiral galaxies). The average fractal dimension in this case is {D^) = 1.81 with a 
standard deviations of a = 0.14, and the average Ohoot is ~ 0.03. By considering the 
associated uncertainties we can say that the values do differ among themselves. 



3.1 Dependence on galactic properties 

To examine possible correlations we propose, in a first approximation, a linear 
model linking all variables and we evaluate its "goodness" via the Akaike's In- 
formation Criterion |1J. This criterion not only rewards goodness of fit, but also 
includes a penalty that is an increasing function of the number of free parameters 
to be estimated (discouraging overfitting). Then, starting from a linear regression 
model based on n variables we apply a stepwise regression process to select the 
subset of variables providing the best model. We begin with the variables T , A, Mb, 
R25, Vrot, and another variable that we have called the average surface density of 
star forming regions (A^(HII)/i?25)- The appHcation of this procedure using the R 
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Fig. 4 Positions of the HII regions in NGC 6946 for high, medium, and low brightness regions, for 
which the resulting fractal dimensions were Dc = 1.64, = 1.82, and = 1.79, respectively. 

environment for statistical computing lfT6l yields that Mg and 7?25 represent the best 
linear model fitting the data: 

D,. = -(0.069 ±0.025)Mb- (0.006 ±0.003)/?25 + 0.456 . (2) 

The overall fit is significant at the confidence level of 95% and the regression coef- 
ficients of the variables Mg and /?25 significant at the confidence levels of 99% 
and 95%, respectively. 

According to our results Dc correlates with Mb and, to a lesser extent, with /?25- 
Brightest galaxies have fractal dimensions higher than faintest ones, while Q find 
a pure scatter diagram when comparing Dc and Mg which, together with the lack of 
correlation with other properties is used as argument favoring a universal Dc value 
in spiral galaxies. In contrast, ifTSl and ifTSl report the same trend in their samples 
of dwarf galaxies, i.e. the highest fractal dimensions for the brightest galaxies. 



3.2 Possible evolutionary effects 

We have analyzed the dependence of Dc on the brightness of the HII regions. We 
first selected galaxies from the sample for which we have data on HII region bright- 
ness and a sample size of A^(HII) > 600. There are 9 galaxies in our sample fulfilling 
these requiriments. We divided them in three equal subsamples ordered in descend- 
ing brightness: the brightest HII regions, the medium bright ones, and the faintest 
ones. Figure |4] shows the resulting distributions in the galaxy NGC 6946. It can be 
seen that the obtained Dc values are consistent with the appearance of the HII region 
distributions. The high brightness map has a relatively small Dc value (~ 1.6) that 
corresponds to a more irregular distribution having clumps and/or filaments sepa- 
rated by low density (or empty) regions. The medium and low brightness maps have 
higher Dc values (~ 1.8) corresponding to more homogeneous distributions. We re- 
calculated Dc for each case and found that for 7 of the 9 galaxies Dc is smaller for 
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the brightest HII regions than for the rest of the data. The brightest regions should re- 
flect, in a first approximation, the initial distribution of star-forming regions in each 
galaxy. It seems likely from our results that some kind of evolutionary effect tends 
to randomize (homogenize) in some degree the initial distributions of HII regions. 



4 Final remarks 

There are many other galactic properties from which Dc could depend, such as star 
formation activity, mass, age, metallicity, or a combination of them. It has been 
suggested that the fractal dimension of the distribution of star-forming sites could 
be increased during the star formation process [4] . Recently, a transition scale from a 
lower to a higher correlation dimension has been suggested as a possible explanation 
for the wide range of observed values 1141 . Obviously, a more complete analysis 
including more galactic variables and a wider and more diverse sample of galaxies 
would be necessary in order to obtain a clearer picture. 
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